Effect of Electric Field on Diffusion in Disordered Materials 
I. One-dimensional Hopping Transport 
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An exact analytical theory is developed for calculating the diffusion coefficient of charge carriers 
in strongly anisotropic disordered solids with one-dimensional hopping transport mode for any de- 
pendence of the hopping rates on space and energy. So far such a theory existed only for calculating 
the carrier mobility. The dependence of the diffusion coefficient on the electric field evidences a 
linear, non-analytic behavior at low fields for all considered models of disorder. The mobility, on 
the contrary, demonstrates a parabolic, analytic field dependence for a random-barrier model, being 
linear, non-analytic for a random energy model. For both models the Einstein relation between the 
diffusion coefficient and mobility is proven to be violated at any finite electric field. The question on 
whether these non-analytic field dependences of the transport coefficients and the concomitant viola- 
tion of the Einstein's formula are due to the dimensionality of space or due to the considered models 
of disorder is resolved in the following paper [Nenashev et al., arXiv:0912.3169 , where analytical 
calculations and computer simulations are carried out for two- and three-dimensional systems. 

PACS numbers: 72.20.Ht, 72.20.Ee, 72.80.Ng, 72.80.Le 
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I. INTRODUCTION 

Charge carrier transport in disordered materials - in- 
organic, organic and biological systems - has been in the 
focus of intensive experimental and theoretical study for 
several decades due to various current and potential ap- 
plications of such materials in modern electronic devices 
(see, for instance, Ref. |l| and references therein). An 
essential part of the research is dedicated to studying 
the mobility of the charge carriers, [a, and their diffusion 
coefficient, D, as the decisive transport coefficients re- 
sponsible for the performance of most devices. Among 
other features, the relation between these two trans- 
port coefficients is the subject of intensive research, since 
this relation (called the "Einstein relation") often pro- 
vides significant information on the underlying transport 
mechanismi 1 - In numerous experimental studies on or- 
ganic disordered materials, essential deviations from the 
conventional form 



(1) 



of this relation have been recognizedi^i&i^&L^ In Eq. (fT]) 
e is the elementary charge, T is temperature and k 
is the Boltzmann constant. Einstein derived this re- 
lation between fi and D for the case of thermal equi- 
librium in a non-degenerate system of charge carriers. 
Deviations from Eq. |T]) were predicted theoretically 
for non-equilibrium transport at low temperature o 9 ' 10 ! 11 
and also for equilibrium transport in degenerate sys- 
tems if the density of states (DOS), which can be 
used by the charge carriers, strongly depends on en- 
ergy, for instance, exponentially 12 or according to a 



Gaussian distribution.— Usually the former DOS is as- 
sumed for inorganic amorphous semiconductors, while 
the latter one is assumed for disordered organic ma- 
terials, such as molecularly doped and conjugated 
polymers i 14 d 5 i 16 i 17 i 18 i 19 In this paper and in the following 
one (Ref. H3) we derive general equations for calculating 
the diffusion coefficient and the mobility of charge carri- 
ers and apply them to systems with the Gaussian DOS, 
since most of the experimental evidence for the violation 
of Eq. |T]) has been reported for organic disordered ma- 
terials. The DOS is taken as 
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where N is the spatial concentration of conducting states 
and a is the energy scale of the DOS distribution. 

Remarkably, experiments on disordered organic ma- 
terials evidence that at relatively low electric fields, at 
which the carrier mobility \i is field-independent and 
hence the carrier transport can be treated as Ohmic one 
(low-field regime), the diffusion coefficient D of charge 
carriers and concomitantly the relation between \x and 
D become essentially dependent on the magnitude of the 
applied electric field F£££& Our aim in this paper and in 
the following one^ is to provide an analytical theory for 
the field-dependent diffusion coefficient and mobility of 
charge carriers. The theory will be checked by computer 
simulations. 

Charge transport in disordered organic materials is 
dominated by incoherent hopping of electrons and holes 
via localized states randomly distributed in space, with 
the DOS described by Eq. © ( J4,lgJ&iL18 1 l& The tran- 
sition rate between an occupied state i and an empty 



2 



state j, separated by the distance r^, is described by the 
Miller- Abrahams expression^! 

Tij ;=voe- 2 ^ I ' >° , (3) 

3 \ 1 , Ae y < y ' 

where vq is the attempt-to-escape frequency. The energy 
difference between the sites is 



£ ij — Ej £ i C X j X i 



(4) 



where the electric field F is assumed to be directed along 
the ai-direction. The localization length of the charge 
carriers in the states contributing to hopping transport 
is a. We assume the latter quantity to be independent 
of energy and we will neglect correlations between the 
energies of the localized states, following the Gaussian- 
disorder-model of Bassler i 17 i 18 i 19 

The challenging problem arises of how to describe the- 
oretically the field-dependent diffusion of charge carriers 
in the hopping regime within the Gaussian DOS. This 
very problem was addressed in the numerical simulations 
by Richert et al. 22 Using a Monte Carlo algorithm with 
a randomly distributed parameter a (the so-called off- 
diagonal disorder), it was shown that the diffusion co- 
efficient for hopping transport in the Gaussian DOS de- 
pends essentially on the field strength at such low electric 
fields that the mobility of charge carriers remains field- 
independent^. This result was interpreted in analytical 
calculations by Bouchaud and Georges,— who considered 
a hopping process in a one-dimensional (ID) system of 
equidistant localized states with transition rates essen- 
tially different from those given by Eq. ([3]) . In the calcu- 
lations of Bouchaud and Georges^ 3 - the transition rates 
between the neighboring sites were taken as 



T v ,±i = r exp 



Ai±i,j ± eFd 



2kT 



(5) 



with Ai t i + i = A J+ i,i distributed according to g(Aij) 
given by Eq. |2]). We will call this model the Random- 
Barrier-Model (RBM) in contrast to the model described 
by Eqs. (J2J) and ©, which we call the Random-Energy- 
Model (REM). Bouchaud and Georges 2 ^ suggested for 
the field-dependent part of the diffusion coefficient in the 
RBM the expression D(F) - D(0) oc Fexp[3cr 2 /8(fcT) 2 ], 
which they claimed to be precisely the dependence found 
in Ref. |22|. Later the authors of Ref. [22] studied the 
quantity D(F) — D(0) by computer simulations in more 
detail 2 ^ and found a quadratic dependence of D(F) — 
D(0) on F at low fields and no turn-over to a linear field 
dependence as suggested by Bouchaud and Georges. 23 
The question arises then on whether this discrepancy in 
the field dependences of the diffusion constant between 
the computer simulations 2 ^ and analytical calculations^ 3 - 
is due to different models (RBM 2 ^ against REM 2 ^), or 
it is due to different dimensionalities considered in these 
two approaches (ID in analytical calculations 23 against 
3D in computer simulations 24 ). The only way to answer 



this question is to obtain exact results for the REM in 
f D and to compare them with the results for the RBM in 
ID on one hand and with the results for the REM in 3D 
on the other hand. This task demands developing a new 
analytical method for calculating drift and diffusion in 
ID systems for the hopping transport mode. In Sec. |TT] 
we present such method. We also present in Sec. IIIII 
the exact result for the field-dependent diffusion in the 
RBM, which differs from the one given by Bouchaud and 
Georges^ Sec.|TV]is devoted to analytical results on the 
field-dependent diffusion coefficient and mobility in the 
REM in the fD case. The exact results for both RBM 
and REM give a linear field dependence of the diffusion 
coefficient at low fields. In Sec. [V]we present the results 
obtained by computer simulations in ID systems. Con- 
cluding remarks are gathered in Sec. I VII 

The following paper— is devoted to diffusion in 3D 
systems. The results in the 3D case clearly demonstrate 
a quadratic field dependence of the diffusion coefficient 
at low fields. One should then conclude that the dis- 
crepancy between the linear— and the quadratic* 2 ^ field 
dependences of the diffusion constant reported in the lit- 
erature is due to the different space dimensionalities con- 
sidered in the two approaches. One should note that 
the differences between ID systems and 3D systems with 
respect to the field-dependent diffusion coefficient have 
been reported in the literature, albeit for systems with 
essential correlations between energies and spatial posi- 
tions of localized states involved into the hopping trans- 
port. Relying essentially on such correlations, Parris et 
al£^ obtained an exact result for the field-dependent dif- 
fusion coefficient in ID systems, which was not confirmed 
in computer simulations carried out on 3D correlated 
systems^, Our study leads to a similar conclusion for 
the Gaussian disorder model without space-energy corre- 
lations. This study was necessary, since the theory from 
Ref. cannot be applied to the case of uncorrelated 
disorder. 



II. ANALYTICAL METHOD 

This section is devoted to one-dimensional hopping in 
the presence of an electric field. The considered system 
consists of a chain of sites separated by a constant dis- 
tance d. Each site is either empty or occupied by a car- 
rier. We consider the limit of small carrier concentration, 
therefore the probability for the ith site to be occupied, 
Pi, is small for each i. The time evolution of probabilities 
Pi is described by equation 

Qp^ 

tt~ = ^i-i,iPi-i +T i+ i ,iPi+i - (Ti i^ 1 +Ti i +1 )p i , (6) 
at 

where is the rate of transition from site i to site j. 
Transition rates are assumed to be time- independent; 
to be non-zero only for nearest neighbors (Tij ^ <^> 
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— — 1); and to obey the principle of detailed balance: 

(7) 



r M+ i £j - e l+ i + eFd 

— exp — 



kT 



where Si is the energy of a carrier on the ith site without 
the electric field, and F is the electric field strength. 

Our aim is to obtain analytical expressions for diffusion 
coefficients with transition rates Tij chosen according to 
either RBM or REM. A similar problem was considered 
by Derrida— who obtained exact results for diffusion co- 
efficient in finite systems with arbitrarily chosen transi- 
tion rates. But, in the limit of an infinite system, his 
expression (Eq. (47) of Ref. [27| ) contains an uncertainty 
of type "oo — oo" , and resolving this uncertainty is a non- 
trivial task. Derrida considered an infinite system only 
for the case if T^ are random independent variables, ex- 
cept that only Tij and Tji may be correlated. This con- 
dition is fulfilled for the RBM, but not for the REM, in 
which Tij and Tjk are correlated due to the common site 
j. Therefore Derrida's approach can hardly be general- 
ized to for the REM. Here we propose another analytical 
approach for evaluating the diffusion coefficient in the 
infinite disordered one-dimensional systems. Derrida's 
method uses a definition of the diffusion coefficient D 
related to random walks: 



- lim — 

2 i^oo dt 



D ~ lim — ({x 2 {t)) 



(8) 



where x(t) is the position of the particle at time t. On 
the contrary, our method is based on the macroscopic 
definition of D as a ratio of current flow j and the long- 
scale gradient of the concentration n of particles: 



D 



dn(x)/dx 



(9) 



We believe that both methods give the same results, 
though our method has an advantage of providing an 
explicit expression for D in the general case of the in- 
finite one-dimensional system (see Eqs. (|27p. (|2"9")l . and 
([55)) below). This expression can be straightforwardly 
applied to the particular cases of the RBM and REM. 

We start by considering the continuous-medium ap- 
proximation. This approximation deals with the carrier 
concentration n(x,t) averaged upon a sufficiently large 
spatial scale. The time evolution of this concentration 
obeys the Fokker-Planck equation 



dn 
~dt 



dn 
dx 



D 



d 2 n 
dx 2 ' 



(10) 



provided that n varies in space sufficiently slowly (that 
is, the characteristic scale of spatial variation is large as 
compared to the scale of averaging). Here v is the drift 
velocity and D is the diffusion coefficient. Let us consider 
the initial concentration n{x 1 0) in the form 



with an infinitely small factor rj. The solution of Eq. (fTO]) 
with the initial condition (fTTj) reads 



where 



i(x, t) — hq exp(?7a; — \t), 



A = vrj — Drj 2 . 



(12) 



(13) 



Since both A and rj are infinitely small, one can resolve 
Eq. (|13|) with respect to rj in the following way: 



r, = -A + -^A 2 + 0(A 3 



(14) 



We will use Eq. (fl"4"|) for calculating the drift velocity v 
and the diffusion coefficient D. For this aim, we need 
a microscopic definition of the coefficients A and rj ex- 
pressed in terms of occupation probabilities pi rather 
than in terms of the concentration n. 

To obtain an exponential time dependence of the con- 
centration, n ~ exp(— At), we can simply postulate that 
each probability pi depends on time in the same way, 
Pi ~ exp(— At). Therefore dpi/dt = —Xpi, and Eq. (J6j) 
can be written as 

-Xpi = Ti-i ti pi^\+Ti + ijPi + i — (ri t i-i+Ti,i + i)pi. (15) 

This is the way of introducing A on a microscopic scale. 

For the spatial dependence of probabilities, one cannot 
expect an analogous form, pi ~ exp(r)di), if the system 
has spatial disorder, i. e. no translation symmetry. In- 
stead, we expect that 



Pi = poQexp^di), 



(16) 



where the coefficients Ci does not exponentially grow or 
decay when i tends to infinity. Consequently, 



which gives 



or, equivalently, 



log — =ndi + 0(l), 
Po 



r 1 l Pi 
V = . I™ —log—, 

i— >±oo dl po 



i A Pi+i 

a \ Pi 



(17) 



(18) 



(19) 



where angle brackets denote averaging over the site num- 
ber i. 

Eq. (|19|) can serve as the microscopic definition of r\. 
However, it is more convenient for our aim to define r\ in 
another way: 



1 /, JM+i 
r) = - { log - 



(20) 



n(x, 0) = no exp(nx) 



(11) 



d \ Ji-i 

where is the flow of carriers from site i to site i + 1 

= Tm+i Pi - Tj+i/i Pi+i- (21 
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It is easy to show that Eqs. (fl9|) and (|20[) give equal 
values of r\. Indeed, in a macroscopic consideration the 
flow of particles j(x, t) is connected to the concentration 
n(x, t) as 



j = vn — D dn/dx. 



(22) 



Therefore, if n ~ exp(7ya:) then j ~ exp(Tyx). Going to 
a microscopic picture, one can get Eq. (|19|) from n ~ 
exp^x) and Eq. (|20|) from j ~ exp^x). Consequently 
the value of r\ should be the same in all these equations. 

Let us now obtain v and D from Eq. (|20[). For this 
purpose we rewrite Eq. (fT5)) taking into account Eq. (f2"T) : 



Aw = ii-i,t - j 



which gives 



Ji.i+1 ^ Pi 

3i — l,i Ji— 1,2 



(23) 



(24) 



The ratio Pi/ji—i,i is a function of A since the probability 
Pi and the carrier flow j are defined by a A-dependent 
equation (115p . We expand this ratio in a Taylor series: 



Pi 



Ji-l,i 



(H + Xb t + 0{\ 2 ). 



(25) 



(Our coefficients a,i are the same as Derrida's r n in 
Ref. H3.) Substitution of Eqs. (J2U) and IJ25J into Eq. (PZOD 
gives: 



^i/logjl + A/i- 



i(log(l + Aa, + A 2 fe J +0(A 3 ))) = (26) 
^(a l ) + ^((fo l )-( a ?)/2)+0(A 3 ). 



Comparing the latter equation with Eq. (|T4]) . one ob- 
tains v and D: 



«=7^, D = * {bi) 7 { $ /2 - (*) 



The expression for v coincides with that obtained by Der- 
rida (Eq. (63) of Ref. l27t ). whereas the expression for D 
is a new result. 

In the rest of this section, we obtain explicit expres- 
sions for the quantities and bi. The mean values (at), 
(a 2 ), and {bi) will be evaluated in Sec. [Hi] for the RBM 
and in Sec. HVI for the REM leading to the analytical ex- 
pressions for the diffusion coefficient D in the RBM and 
in the REM. 

In order to find the coefficients a,, we set A to zero 
in Eq. ([25]) . As it is seen from Eq. (j23|) . the carrier flow 
does not depend on i in the case of A = 0. Dividing 
Eq. (|21[) by the carrier flow, one obtains a set of equations 
for coefficients a,: 



The solution of Eq. ([28)1 can be presented as an infinite 
series: 



1 



ri+i,irj + 2 i i+i 



+2,i+3 



(29) 



what can be checked directly by substituting Eq. (|29|) 
into Eq. (|28p . To prove the convergence of the series (|29p , 
let us rewrite it using the condition of detailed balance, 
Eq. 0: 



a t = Trj +1 + B- 1 exp (^i) rr + \ 

B-^f^r^i..., 



i+2 



(30) 



where B = exjp{eFd/kT). For any physically reasonable 
system, the quantities exp ((s i+k - Zij/kT) T~^ k i+k+1 
can be regarded as having an upper boundary. Denoting 
this boundary as C, we get an upper estimate for af. 



a t <C + B~ 1 C + B~ 2 C+-' 



C 



1 - B' 



(31) 



that proves convergence of the series (|2"9")l under the con- 
dition B > I, i. e., eF > 0. 

In order to obtain bi, we need a set of equations con- 
necting bi to in analogy with Eq. (|28|) that connects 
di to a^+i. We will derive the necessary equations using 
Eq. ([2"TJ). Eq. (|2"3]1. and the Taylor expansion Let 
us first divide Eq. (|2~lj) by j'^i+i and slightly rearrange it: 



p Pi r Pi+i _ -, 

-L : -L i+l.i- — J-- 



(32) 



Ji,i+1 Ji— l.i Jiji+l 

Let us now use the expansion (|25|) for quantities Pij 

(Oj + Abj) - r i+ i^(ai + i + A6j + i) = 1 + C(A 2 ). 

Ji,t+1 

(33) 

The latter equation contains the ratio ji—ij/jii+i- We 
derive this ratio from Eq. (|24| using also the expan- 
sion 



3i-l, 



I 



1 



1-A ai +0(A 2 ). 
(34) 



1 + Xpi/ji-i.i 1 + Aa, + 0(A 2 ) 
Finally, let us substitute Eq. (fM)) into Eq. 



r M+ i(l-Aa l )(a J +A&. i )-r 4+ i.. i (a. 1+ i+A& i+ i) = 1+0(A 2 ), 

(35) 

and collect separately terms, which do not contain A, 
and those proportional to A. The former terms lead to 
Eq. (|2"8|) , while the latter ones give the equation 



T iti+1 (Xbi - Xaf) - Y l+ i.i\b i+ i = 0. 



(36) 



Eq. (f3"6"| provides a desired set of equations for coefficients 
bi-. 



Vi r^i+i a, — IYfi,^ ai+i — 1. 



(28) 



v« r^i + i 6^ — rj_|_i^ 6j_|_i — r^i + i a 2 . 



(37) 
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The solution of Eq. (j3T[) can be found as an infinite 
series: 



hi = a, ; 



2 



-a 



i+M 1 i+2,i+l 2 



i+1 



i,i+l-L i+l,t+2 



*i+2 



(38) 



that can be checked by substitution into Eq. ([57)1 . Like 
Eq. ([29"|). the series l[3"8"|) converges provided the product 
ei* 1 is positive. To prove it, we substitute the condition 
of detailed balance, Eq. ([TJ, into this series: 



h = 
B~ 



' exp 



- B 1 exp 

/ e i+2 -Si \ 

\ kT J 



( 1+ fcT t ) a i+l" 



(39) 



'i+2 



In any real system we find an upper limit for the quan- 
tities exp ((e i+k - Si)/kT) 

a i+k~ Setting this limit equal 
to C, we obtain an upper estimate for 6^: 



h <C + B- 1 C + B- 2 C + --- = 



C 



1 - B- 



(40) 



Therefore the series (|38|) converges if B > 1, i. e. if 
eF > 0. 

As a result, we have obtained an analytical ex- 
pression (|27|) for the diffusion coefficient D in a one- 
dimensional hopping system. For coefficients a% and &i 
that contribute into Eq. (|27[) we have found series repre- 
sentations (f29|) and (l38|) in the case eF > 0. It is easy to 
write down analogous series for and hi in the opposite 
case, eF < 0. 



RANDOM-BARRIER MODEL: EXACT 
RESULTS 



III. 



Let us now apply Eqs. (J27J, (|29"]1. and d35J) to the 
random-barrier model described by Eqs. ((2|) and ([5]). In 
this model, any two transition rates and are sta- 
tistically independent, if (ij) and (kl) are different pairs 
of sites. The rates Ti^+i and T^i^, related to the same 
pair are connected to each other. As a result, all statisti- 
cal properties of the random-barrier model are defined by 
mean values (r™ +i r" +1 J with different m's and n's. Wc 
introduce the following notations for these mean values: 

"ii = (rj+i.j/r^i+i), 

™ 2 = (rf +M /r? <+1 >, 

m 3 = (l/T iti+1 ), (41) 

m 4 - (i/r? i+1 ), 
m 5 = (r i+M /r? i+1 ). 

In order to obtain the drift velocity w and the diffu- 
sion coefficient D from Eq. (|27p , one should calculate the 
mean values (a^), (a 2 ), and We start with calculat- 
ing (cii). Let us denote successive terms in the expan- 



sion J29l) as a (0) , a (1) , a (2 ), 
1 



Then, 



,(°)\ = 



,(2) 



M,t+1 
M,*+l 
M,i+1 



= m 3 , 



1 



M+2,t-|-l 



+ l,z+2 



= mim 3 , 
1 



(42) 



mfm3, 



k 



Consequently, 

(a,:) = (a (0) ) 



,(2)\ 



7713(1 + 777-1 + mf + . . .) = 7713/(1 — 77ll). 



(43) 



The mean value (a 2 ) can be represented as a sum of 
values (a^a^) over all pairs k, I: 



00 00 



EE<« (fc) « 

fc=0 1=0 



(44) 



It is easy to check that 



i\ k 1 m2 , m3m5, if k < I, 
ifk = l, 



77737774 



Then, presenting Eq. (|44|) in the form 



(45) 



' 1 m2m 3 m 5 , if k > Z. 



Wad) 



K 2 > = E« a(fe) ) 2 ) + 2 E< a 



(46) 



fc<; 



using Eq. (|45[) and introducing the notation p = Z — k — 1, 
we obtain 

OO CO 00 

(of) = J] 777^7774 + 2 J] J] m[~ k ^ 1 m^m 3 m 5 = 

k=Q k=Ql=k+l 



— + 2777 3 m 5 V 777$ V 777? = , . 

1-IB2 13 ^ ^ ^ 1 1— 1T12 (1— mi 1- 1712) 

fe=0 p=0 

(47) 

In an analogous way, we denote successive terms of the 
series (|3"8")l as b^-°\ b^ 2 \ . . . The mean values of these 
quantities are 



2to 3 to 5 



(b {0) ) 
(& (2) > 



/r. i+ i,i 2 \ 2 

( y~ -o-i+i ) =mi{a i ), 



Tj+l. ^1+2.1+1 2 

a. 



( + 2 



mi{aj), (48) 



(&«)=m*(a 2 ). 
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Therefore, 

(6i) = <6( )) + (6Ci)) + (6{2)) + ... = 

(a?) (1 + m x + m? + . . .) = (o?)/(l - mi). 



(49) 



Finally, we substitute Eqs. (|43|). (|47j) . and (J49J) for the 
mean values (a*), (a?) and (bj) into Eq. (|2"T)) . This gives 



= <*■ 



1 — mi 

W3 



(50) 



D 



d 2 



2mJ(l 



m 2 J 



(7714(1 — mi) + 2TO3TO5). (51) 



Equations ([50)1 and (|?T|) are not new results — they were 
obtained in Ref. [2JJ (Eq. (67) and Eq. (70), respectively). 
Their derivation in the frame of our method clearly 
demonstrates that both methods (Derrida's and ours) are 
consistent. In the rest of this subsection, we will apply 
these equations to the case of Gaussian distribution of 
barrier heights. 

For the transition rates defined by Eqs. ^ and (J5J), 
the mean values mi . . . 777,5 are easy to evaluate. Setting 
Fq equal to unity for the sake of simplicity, we obtain 



mi 


= exp ( 


' eFd\ 
, kT ) ' 


m 2 


= exp ( 


' 2eFd\ 
kT ) ' 


m 3 


= exp | 






m,4 


= exp ( 


f 


K 2(kT)^ 


m 5 


= exp | 


' 0* 


1 8(kT) 2 



(52) 



eFd \ 
2kT J 

eFd 
kT 

3eFd 
2kT 

From Eq. (f5"0"| one obtains the result for the drift velocity 
v: 

- 2 ^ 'eFd s 



v 



2d exp 



sinh 



2kT 



(53) 



Note that we derived the latter equation only for the case 
eF > 0. However, it is easy to show that Eq. (|53^) is valid 
for any direction of the electric field. Indeed, the right- 
hand side of the equation is an odd function of the electric 
field F. The left-hand side (drift velocity) should also be 
odd, because the system is symmetrical with respect to 
a left-to-right mirror reflection [x — ► —x, F — ► — F, v — > 
—v). Therefore, if Eq. (|53[) is satisfied for positive electric 
fields, it remains valid for negative fields, and vice versa. 

An expression for the diffusion coefficient D as a func- 
tion of F can be obtained by substituting the mean val- 
ues (|52|) into Eq. (|5Tj) . Strictly speaking this procedure 
is valid for D only in the case eF > 0. One can how- 
ever generalize this expression for any sign of the electric 
field using the fact that (for symmetry reasons) D is an 
even function of F . One simply should replace eF by its 
absolute value, \eF\ in all expressions. The result reads: 

r, ,2 ( c 2 \eF\d\ 

D = d Hp [-w"w) + 

\8(kT) 2 J \ 2kT J 



(54) 



Eq. (|54p was obtained for non-zero electric fields. How- 
ever, one can check that it holds also for F = 0. 



Eq. (|54|) differs from the expression given by Bouchaud 
and Georges A D(F) - D(0) cx Fexp[3cr 2 /8(fcT) 2 ], 
though it is linear in F to first order. Eq. (154")) is plot- 
ted in Fig. [21 together with numerical results obtained in 
SecM 



IV. RANDOM-ENERGY MODEL: EXACT 
RESULTS 



The random-energy model in one dimension implies 
the following definition of transition rates: 



i,i±\ 



exp 



Ae M± i + |Ae M± i 
2kT 



(55) 



where Ae^i-ti = Si±\ — Si =F eFd is the difference between 
the energies of a charge carrier on the final site and on 
the initial site, respectively for each jump. For simplicity, 
we set the constant To = vq exp(— 2d /a) to unity. 

For the REM, one can use the same way of calculating 
the velocity and the diffusion constant as for the RBM. 
The REM contains more correlations between transition 
rates that the RBM, which leads to more complicated 
calculations of the mean values (cij), (a 2 ), and (fej). In 
the REM, each rate depends on the energies £{ and 
Ej, which are independent random variables. Therefore, 
the rates and Tki are correlated if the pairs of sites 
(ij) and (kl) have at least one site in common. 

We will see below that the drift velocity v and the diffu- 
sion coefficient D depend on eleven quantities mi . . . mn 
related to the statistics of site energies and transition 
rates: 



mi = 




m 2 = 


( e -2 6 ./*T) j 


m 3 = 




ni4 — 




m 5 = 


(e- ei/fcT r^ l); 


m 6 = 




7717 = 




m 8 = 




m 9 = 


/ e (E<+l-E<)/*r r -) 


mio = 


(e^/^T-^T, 


mu = 


^ e (e i+1 +Ei)/kT p-] 



(56) 



i+2/> 



+ 1 ^ i+1,1+2/- 



In the following we proceed for the REM along the 
same steps as for the RBM in the previous section. 
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1. Calculation of (at) 
Let us denote successive terms of the expansion ([30 



asa<°\ oW, a( 2 ), and so on. Then (a,) = (a^) + (a^) + 
(a^ 2 )) + . . . . The latter quantities can be easily expressed 
via mi, m3 and 7774: 

<a (0) > = <r M Vi>=™3, (57) 

( «) = B-V^'^w) = 

(fl ( 2)) = ^/^^/B.^ = B - 2miTO4j 

and, generally, (a^ 5 )) = B~ k m\m± for any fc > 0. 
Then, 

(04) = m 3 +( J B" 1 +S- 2 +. . ,)m 1 m 4 = m 3 + ^— 7. (58) 

75 — 1 



U. Calculation of (af } 

According to Eq. ([44|) . the calculation of (af) is re- 
duced calculating the mean values (a^a^) for all inte- 
ger k > and Z > 0. Thus one can reduce the mean 
values (a^a^) to: 



(a^a^) 



if k = 1 = 0, 
if k = 0, I = 1, 
if fc = 0, Z > 1, 
if fc > 0, I = k, 
B- 2k - 1 m 2 m 11 , if k > 0, Z = k + 1, 
if fc > 0, Z > fc + 1. 



m 6 , 

B~ l m4m5, 
B~ 2k mim%, 



B k l m2m\ 



(59) 

The next step is the estimate of the infinite series (14"4"1) . 
It is convenient to rearrange the summation in Eq. (|44[) , 
separating terms corresponding to different lines of 
Eq. (59): 



(60) 



(of) = ((a(°)) 2 ) + £ ((a( fe )) 2 ) + 

fc=l 

OO 

2(tt(°)a( 1 ))+2E(a(°»a( | ))+ 

i=2 

OO 00 00 

2 £ (aW a ( fc+1 )) E (a (fc) a( fc +f)), 

fe=l fc=lp=2 



where p = I — k. Substituting Eq. (|59[) into this expan- 
sion, one gets: 



(a 2 ) = m 6 + m 2 m 8 B 2fc + 2B 1 m 9 - 

k=l 



21714,1715 B 1 + 2to2«7ii J]] B 



2k-l 



(61) 



i=2 



fc=l 



2m 2 mlY J B- 2k Y J B ~ P - 

k=l p=2 



Finally, one can sum up the geometric series. The 
result reads: 

_,2\ m _ 1 2m;j 1 2yn4m.5 , 



m 2 



m 6 



1 B(B-l) 1 

2m2rri4 
B(B-1)(B 2 -1) 



(62) 



5. Calculation of (bi) 

In an analogous way, (bi) can be expressed as a sum 
(6(0)) + ( & (i)) + (ft(2)) + 5 where 6 (0) j 6(1)... 6(0 } ... are 

the terms of the expansion ([39]) . Keeping in mind that, 
according to Eq. (f2"9"| . the values aj+i, ■ ■ ■ do not 

depend on £j, one can express the mean values M fe ) as 
follows: 

(&W>) 

(6 (1) > 
(6 (2) > 



(a 
t3 



1/ 



B- 2 (e 



-Ei/fcT 
-Ei/fcT 



/kTJ 



— /^a 2 +2 ) 



and so on for larger i, where M 



= B^rmM, (63) 
= B- 2 mi M, 

(64) 

(e^/ feT a 2 ). Thus, 



(6i) 



)miM, 



or 



(6*) - (o?) + 



mi 
B — 1 



■M. 



(65) 
(66) 



In order to find M, one can expand it in series analo- 
gous to Eq. dH}: 



(67) 



k=o z=o 

and express each term of the expansion via: 



(e&a {k) a {l) ) = < 



7777 




if fc = Z = 0, 




73- 




if fc = 0, z = 


1, 


73- 


'm 4 m 3 , 


if fc = 0, z > 


1, 


73- 


ZK TOlTO 8 , 


if fc > 0, z = 


fc, 


73- 


9h— 1 

mimii, 


if fc > 0, z = 


k + l, 


. 73- 


it— / 9 

TOlTOJ, 


if fc > 0, Z > 


fc + 1. 



(68) 

The following steps are the same as the ones lead- 
ing from Eq. (|59|) to Eq. (|62"jl . Instead of proceeding 
in this way, one can recognize that Eq. (|59| transforms 
to Eq. (f6"8"|) by the following replacements: 



m 6 -> m 7 , 



m 9 



"7io, m 5 -> m 3 , m 2 



7771. 



Applying the same replacements to Eq. (|62|) , we obtain 
the result for M: 



M = m 7 



2mm 
B 



2m^m^ 
B(fl-l)' 



(69) 



m l B(B^-l) + B(B-l)(B^-l) • 

Substituting Eqs. (|6"2")l and (|6"9")l into Eq. (|tj6"l) . one obtains 
the expression for (6j) in terms of the values mi . . . mn. 
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4- Drift velocity and diffusion coefficient 
Combining equations ((27}, ([58]), ([62]), (J66J) , (J69J) leads 



to 



J 



TO.3 + 77117714/ (-B — 1) 



(70) 



D 



2d 



m 6 



2m\ 
B- 1 



7777 



2mg 2777,4777,5 

2mio 
t3 



73(73 - 1) 

277737774 
B(B-l) 4 



777 2 



777i- 



i?777g + 2mn 

t3(t3 2 - 1) 

737778 + 2?77ll 

t3(t3 2 - 1) 4 

I 



277727774 



t3(t3- 1)(t3 2 - 1) 
2?77im| 
t3(t3- 1)(t3 2 - 1) 



(71) 



where B = exp(eFd/kT). Note that we derived these 
equations for the case eF > 0. One can easily generalize 
the equations for the case eF < 0, keeping in mind that 
v is an odd function of F, and D is an even function. 

Eq. ([70)) was obtained previously by Cordes et al. re- 
using Derrida's method,- 27 - while Eq. [71] is a new result. 

Equations ([70]) . ([7T]) are general for the case of the 
REM — their derivation is not restricted by a special 
choice of the density of states or by the choice of the 
relation between the transition rates and the site 
energies Ei, £j. We used only four assumptions: (i) all 
sites are arranged in the line with constant distance d be- 
tween them; (ii) there are only transitions between near- 
est neighbors; (iii) transition rates "forth" and "back" 
(Ti t i+i and r^+i^) obey the principle of detailed balance, 
Eq. ([7]); (iv) energies of different sites are independent 
random variables having the same distribution function. 



Gaussian density of states 



We will now evaluate the quantities mi . . . m-n, assum- 
ing a Gaussian density of states, Eq. ([2]), and the Miller- 
Abrahams transition rates, Eq. ([55]) . This evaluation is 
straightforward, for example: 



777i =Nq 1 fe- £ / kT g(e)de = 

=~ I CX P (~W ~ &) d£ = ex P ( 



'2-KG 

7773 = N^ 2 JfT^g(e 1 )g(e 2 )de 1 de 
1 



27TCT 2 



JJexp 



£2— ei—eFd+\e2— ei—eFd\ 
2kT 



2(feT) 5 



s 2 +e 2 N 

-hy-^ ) de 1 de 2 - 





\ A, 




/ a,b 


b 


A 60 ° 


2 


1 
1 
1 




a 




2V3 



FIG. 1: Integration area for the definition of the function 
F{a,b). 



The latter integral can be evaluated by substitution u = 
£\ + £2, v = £\ — £2, and the result is 



7773 = -erfc 



eFd 
1 / a 2 



eFd 
~~kT 



rfc 



eFd 
~2cT 



a 

kf 



(72) 



where erfc is the complementary error function, erfc (a;) = 
27T- 1 / 2 f™e~ t2 dt= l-erf(z). 

The values 7779 . . . mxi are triple integrals. They cannot 
be expressed in elementary functions, but they can be 
reduced (by a substitution u = e\ + e 2 + £3, V — e\ — £2, 
w = £1 + £ 2 — 2£3) to a function J- (a, b) defined as 



T{a, b) 



1 



e -( x2 +y 2 Uxdy, 



(73) 



where the area of integration A a ^ is shown in Fig. [T] 

The results are collected in Table HI These results ob- 
tained for Miller- Abrahams hopping rates, Eq. (|55[) look 
rather complicated for analytical estimates. Below in 
Sec. [V] we use these expressions from Table [J for numer- 
ical calculations and present the results for the Miller- 
Abrahams hopping rates, Eq. (|55|) . Here we will proceed 
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TABLE I: The values of mi . . . mn for the random-energy model with Gaussian density of states. "MA" refers to using Miller- 
Abrahams hopping rates, Eq. (|55p. "modified MA" — to hopping rates defined by Eq. (|74jl . Other notations: A = exp(<j 2 /(fcT) 2 ), 
B = exp(eFdfkT), a = a/kT, (3 — eFd/a; erfc is the complementary error function; the function T(a, b) is defined by Eq. ((73}. 



Notation 


Definition 


Value (MA) 


Value (mourned MA) 


nix 




T— 

\[A 


VA 


m,2 




A 2 


A 2 


m 3 


<r^i) 


lEerfcC-D+AB-^erfcd -a)] 


1+AB- 1 


m.4, 




l>/Z[erfc(-f - f ) + B _1 erfc(| - f )] 


VAil + B- 1 ) 


m 5 




i^[erfc(-| + f ) +yl 2 B- 1 erfc(§ - ^)] 


VA(1 + A 2 B~ 1 ) 


rn 6 




i [erfc(-f) + A 4 B- 2 erfc(f - 2a)] 


1 + 2AB- 1 +A 4 B-' 2 


7717 




|>/A[erfc(-f - f ) + A 2 £r 2 erfc(§ - ^)] 


VA(l + 2B- 1 + A 2 B-' 2 ) 


m$ 




±,4 2 [erfc(-§ - a) + B- 2 crfc(§ - a)] 


A 2 {l + B- 2 ) + 2AB~ 1 


mg 


\ L - 1 - i,i+l - 1 - i+l,i+2/ 


AF(-3/3, 2a - /?) + A 3 B-' J P(3/3 - 3a, /3 - 3a) 
+^B _1 [ierfc(^) - ^(3a - 3/3, a - /3)] 
+A 4 B- 1 [ierfc(f - 2a) - ^(3/3, /3 - 4a)] 


A(l + A 2 B- 2 + (1 + A^B- 1 ) 


mio 


/pSi+l/ATp-1 p-1 v 
\ c 1 - 1 - i+l,i+2/ 


V^tF(-a-3/3, a- /3) + A i/2 B-' 2 F(3f3-2a, (3-2a) 
+v / 3B~ 1 [|erfc(-f ) - ^(2a - 3/3, -/3)] 
+^ 5 / 2 B- 1 [ierfc(f - 22L)-.F(a + 3/3,/3-3a)] 


-/A(l + AB- 2 + (1 + A^B' 1 ) 




/ p (e»+l+Ei)/fcr r -l p-1 i 
\ e 1 i,i+l 1 i+l,i+2/ 


AF(-2a - 3/3, -/3) + ^4B~ 2 ^(3/3 - a, /3 - a) 
+AB- 1 [±erfc(^£) - T(a - 3/3, -a - /3)] 
+A 2 B^ 1 [ierfc(f - a) - :F(2a + 3/3, /3 - 2a)] 


A(l + B-' 2 + (1 + A)B~ 1 ) 



with analytical calculations based on slightly modified 
expressions for the hopping rates, which allow straight- 
forward analytical estimates. We suggest to use, instead 
of Eq. (|55|) . the following "modified Miller- Abrahams 
rates" : 



r 



i,i±l 



1 + exp 



Ae 



M±i 



kT 



(74) 



where the constant Tq will be set equal to unity for the 



sake of simplicity. The difference between Eq. (|55[) and 
Eq. (fT4"|) becomes negligible when | Ae^ ^-i-! | ^> kT. There- 
fore, for a kT we expect a good agreement between 
results obtained with these two kinds of hopping rates. 



Results for the "modified Miller- Abrahams rates" are 
also shown in Table QJ Substituting them into Eqs. (|70| . 
([TTjl . one can get the explicit expressions for the drift 
velocity and the diffusion coefficient: 30 



J 



B 



2A + B - 1 ' 



(75) 



D 



2 AA A {B - l) 2 + IQA^iB ~ 1) + 16A 2 + 2A(3B + l)(B 2 - 1) + (B - 1) 3 (B + 1) 

2(B + 1){2A + B-1) 3 : 

I 



(76) 



where A = exp(cr 2 / (kT) 2 ) , B = exp(\eF\d/kT). 

Now we will consider the mobility ^(F) = v(F)/F 
and the diffusion coefficient D(F) in the limit of small 
field. We restrict ourselves to the case of modified Miller- 
Abrahams rates. For F = 0, one can obtain from 
Eqs. ([75]) and 176]) : 



D(0) = M (0) ^ = y exp(-(a/fcT) 2 ) 



(77) 



For small temperatures, kT <C <r, the mobility and 
the diffusion coefficient can be approximated by simple 



expressions: 



ed 2 \Fe\ed 3 F 2 e 3 d 4 



D(F) 



2A 



2AkT AA{kT) 2 12A{kTf ' 

u \Fe\d 3 F 2 e 2 d 4 A \Fe\ 3 d 5 A 
L 2kT + 8{kT) 2 + 16(fcT) 3 



(78) 



(79) 



These approximations are valid for sufficiently small 
fields, \eFd\ < kT. 

From the latter approximated expressions it is obvious 
that there is a cusp at F = for both n(F) and D(F). 



10 



— 

■a 



14 
12 
10 

8 
6 
4 
2 
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X* 
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2 
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■B □■ 
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1 

eFd/o 



1.5 



FIG. 2: £>(F) in the RBM for different temperatures T. The 
curves show the analytical solution Eq. (|54|l . while the sym- 
bols show numerical results (Eq. I80|l for chains with N = 10 7 
sites. The inset shows the low-field behavior. 



T3 




eFd/o 

FIG. 3: D{F) in the REM for different temperatures T. The 
curves show the analytical solution (Eq. (|7ip with mi . . . mu 
from the MA column in Tabled . The symbols show numerical 
results fEg. lSOp for chains with N = 10 8 sites. The inset shows 
the low-field behavior. 



The dependence D(F) demonstrates a linear behavior 
for very small fields (|eFd| <C kT/A), when the first two 
terms in Eq. (|T9")) are dominating, and a parabolic behav- 
ior for intermediate fields (kT/A <C \eFd\ <§; kT), when 
the third term is dominating. 



function of all the hopping rates Ti±ij in the chain^ 7 -: 

iV N 

1 n+i~ 



D 



1 



N 

L n,n+l ^"nXn 



n—l 2—1 



N + 2 
v^^, (80) 



N 



T N r 

l^n=l '« 



jV 



+ l,n 



+1 



n,n+l 



r 



n,n+l 



7V-1 i 

»=i j=i 
»=i i=i 



- n-j + l.n-j 

- n-j,n-j+l 



(81) 



(82) 



(83) 



In this section, we follow Ref. [27] and set the distance d 
between sites equal to unity for the sake of simplicity. 

For both the RBM and the REM, the diffusion coeffi- 
cient was obtained for different temperatures and fields, 
by generating several chains with random jump rates (ac- 
cording to the respective model), evaluating Eq. (|80|) for 
each chain and averaging the results. Long chains (10 7 
and 10 8 sites) were needed to obtain a good agreement 
between different realizations of the chains. 

For chains of this length, a direct evaluation of 
Eqs. ((80|) -([83 |) is not practical. Below, the equations are 
rewritten in a form that can be evaluated in O(N) steps, 
using recursion relations. Define 



and h n 



En,n 



+ 1 



F 



+ 1 



and further 



N-l i 

Gn = 1 + e n gn +i 

i=l j=l 
N-l i 

i=l 3 = 1 



(84) 



(85) 



(86) 



so that r„ = G n / r n>n+1 and u„ = H n / r„ in+1 . All 
G n and H n , and thus r n and u n can now be calculated 
efficiently from 



Gn—i — 9nG n — G + 1, 



(87) 



V. NUMERICAL RESULTS 

In order to verify the analytical results obtained above, 
we perform numerical calculations for a one-dimensional 
chain of N hopping sites using the equations that give 
the drift velocity v and the diffusion coefficient D as a 



Hn+i = h n H n — H + 1 (88) 

where G = H = .91.92 ■ ■ ■ 9n = hihi ■ ■ ■ hjq ■ For the first 
term in the brackets in Eq. (|80p. define S n = J^., ir n+ i 
and S = Y^Li r n- Now 

Sn+i =S n -S + Nr n+1 . (89) 



11 



These relations are numerically stable if G < 1, which is 
satisfied if the average drift is to the right (towards larger 
site indices). The diffusion coefficient is now given by 

„ 1 | H n S n , , 7 v^ G n H n \ N + 2 

(90) 

It seems tempting to write Eq. ([57]) in the form 

G„ = (G n -i +G- l)/g n , 

so that all equations could be evaluated starting from 
n = l, but this form is too susceptible to numerical errors 
to be usable in practice. Thus one has to evaluate all G n 
with Eq. Q87p starting from Gn and store them in a table. 
S n and H n do not need to be stored, since they can be 
evaluated while performing the sum in Eq. (|90[) . starting 
from n = 1. 

With this method of evaluation, numerical results for 
the diffusion coefficient were obtained. The results are 
shown in Fig.[2]for the RBM, together with the analytical 
results, and in Fig. [3] for the REM. For both models the 
diffusion coefficient is linear in the electric field (at low 
fields), see the insets in each figure. 

VI. DISCUSSION AND CONCLUSIONS 

Figures [2] and [3] demonstrate a full agreement between 
analytical results, Eqs. (|54[) and (jTT]) , and numerical ones 
based on Eq. (|80|) . This result can be considered as evi- 
dence that the two definitions of the diffusion coefficient 
D, Eq. ([5]) and Eq. ©, give the same quantity for hop- 
ping in one-dimensional disordered systems. The former 
definition expresses D via the variance of particle dis- 
placement during a random walk, while the latter one 
defines D as a ratio between the particle flow and the 
gradient of macroscopic concentration of particles. Al- 
though it seems to be obvious from a physical point of 
view that both definitions should give the same result, 
no formal proof has yet been known. 

Computer simulations evidence that at low tempera- 
tures the diffusion constant experiences significant fluc- 
tuations from one realization to another, even for sys- 
tems containing millions of localized states. The reason 
of such a large fluctuations in ID hopping is essentially 
the same as in 3D hopping — it is the sensitivity of the 
diffusion coefficient to the very rare sites with low ener- 
gies. We will discuss this phenomenon in detail in the 
following paper.— We also note that sample-to-sample 
fluctuations of mobility for all parameters presented in 
Figures [2] and [3J are negligible. 

The most striking property of the ID diffusion is its 
linear dependence on electric field: 

D(F) = D(Q) + a\F\ + 0{F 2 ) , (91) 

where ci / 0. It means that the diffusion coefficient is 
a non- analytic function of electric field. From general 



physical arguments one can hardly expect such a behav- 
ior. Instead, one can expect that, as D(F ) is an even 
function, it can be expanded in a Taylor series with re- 
spect to F 2 : 

D(F) — D(0) + (3F 2 + 0(F 4 ) . (92) 

The physical reason of this non-analyticity is still un- 
clear. As first steps to acquire an understanding of this 
phenomenon, we will try (i) to provide a mathemati- 
cal explanation of it and (ii) to find out which systems 
demonstrate the linear field dependence of the diffusion 
coefficient and which systems lack this behavior. 

From the mathematical point of view, the possibility 
for a non-analytic dependence D(F) can be seen from 
expressions ([50)1 and (f3"9"| for the coefficients at and bi as 
functions of B = exp(eFd/kT). These expressions are 
series that converge at B > 1 and diverge at B < 1. For 
B < 1 (i. e. for negative eF) one can obtain converging 
series, applying a "mirror reflection" transformation (i + 
k -> i - k, B -> B- 1 ) to Eqs. {3D]) and ([39]). Therefore 
di and bi are defined by different series expansions for 
positive and negative values of the field. Moreover, at 
eF > the values a* and bi depend on quantities related 
to sites i, i + 1, i + 2, ... (for example, on r,; + i^ + 2, 
Ej+i, etc.); on the contrary, for eF < these coefficients 
di and bi depend on a different set of sites: i, i — 1, 
i — 2, ... So, it is obvious that, in a disordered system, 
the function a,i(F) for negative F cannot be obtained 
by analytic continuation of this function for positive F, 
and vice versa. The same is true for bi(F). Keeping 
in mind that the mobility fi = v/F and the diffusion 
coefficient D depend on et's and b's via Eq. ([27]) . one can 
conclude that negative-field parts of the functions fi(F), 
D(F) may not be analytic continuations of their positive- 
field parts. Consequently, a non-analyticity of fJ.(F) and 
D(F) at F = is possible. 

It is evident from Figs. [2] and [3] and from Eqs. ([M]) and 
([79"]) that both the RBM and the REM demonstrate a 
linear low-field behavior of D(F), according to Eq. (|91[) . 
Moreover, Eq. ([75]) shows that the mobility fi(F) in the 
REM also contains a linear contribution with respect to 
F, while in the RBM the mobility is a smooth function 
of F, see Eq. (p]) . 

Parris et al 25 ^ 29 considered a model of a ID continuous 
medium with smooth disorder potential, and obtained 
analytical expressions for n(F) and D(F). Although the 
low-field behavior of the mobility and of the diffusion 
coefficient were not discussed in detail ) 25 ' 29 one can learn 
from equations (25) and (69) of Ref.]29jthat this behavior 
is qualitatively the same as in our REM. The method of 
Refs. [25II29I is, however, not directly applicable to the 
Gaussian disorder model considered here. Therefore a 
separate derivation was necessary. 

One-dimensional transport with Gaussian density of 
states (DOS) has another peculiarity — namely, there are 
sites with arbitrarily high energies, which represent bar- 
riers for transport. Despite the fact that such barriers 
are very rare, their influence on the transport proper- 
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ties should be noticeable (unlike in the 2D and 3D cases, 
where a carrier can easily pass around these hard places). 
One can argue that this peculiarity might be responsi- 
ble for the unusual non-analytical behavior of n(F) and 
D(F). In order to check this assumption, we consider a 
ID system with a discrete DOS g(e) allowing only two 
values of energy (e — — a and e = a): 

g (s) = ls(e + a) + ^S(e-a). (93) 

It is obvious that there are no high barriers in this system. 
Using our general expressions ([50)1 . ([S"Tj) . (fTCTj) . and (fTTj) . 
we have calculated the dependencies n(F) and D{F) for 
this DOS for both RBM and REM. The results demon- 
strate qualitatively the same low-field behavior as in the 
case of Gaussian DOS. Therefore the property of the non- 
analyticity of /i and D cannot be attributed just to the 
presence of infinitely high barriers provided by the cor- 
responding DOS. 

We have also checked whether this non-analyticity is 
related to the assumption of the nearest-neighbor hop- 
ping. We have considered the same random-energy model 
as discussed in Sec. IIV1 except that we allow transitions 
to distant states. The dependence of transition rates 
Tij on distances r^ between the sites is governed by the 
Miller- Abrahams expression ([3]). The mobility and dif- 
fusion coefficient as functions of electric field were calcu- 
lated by a Monte-Carlo algorithm described in the follow- 
ing paper Again, the results have shown linear low-field 
dependencies n(F) and D(F). 

Therefore one can conclude that the linear behavior of 
the diffusion coefficient, Eq. (|9ip. is a robust property of 
ID disordered systems. 

On the contrary, at higher dimensions there are no ev- 
idences of a linear field dependence of the diffusion coeffi- 
cient. Both our Monte-Carlo simulations (see the follow- 
ing paper, Ref. HO) and previous studies 2 ^ clearly demon- 
strate a parabolic field dependence, Eq. (f92|) . for 2D and 
3D systems with Gaussian disorder. It is worth to note 
that a smooth dependence D(F) should be obtained also 
in ID systems without disorder — namely, in systems with 
periodically repeated site energies (and barrier heights). 
Indeed, such a system is equivalent to a finite chain (with 
periodical boundary conditions) considered by Dcrrida. 27 
One can therefore use Derrida's formula for D, Eq. ([80)) . 
This formula, in the case of finite chain length TV, is an 
analytic function of the transition rates Ty, and conse- 
quently of the electric field. Our numerical studies also 
confirm that for small N there is a region of parabolic 
dependence D(F) around F = 0. Thus, disorder is im- 
portant for the non-analytic behavior of D(F). 

Finally, we will discuss the applicability of Einstein's 
relation |T]) for finite electric fields. Our analytical re- 
sults, Eqs. dSS]), ((Ml), ([75]), and (7jj]) show that Einstein's 
relation is violated at any non-zero field F both in RBM 
and REM, and that the deviation from Eq. ([T]) is pro- 
portional to One should note that this phenomenon 
is related to the discrete nature of the systems, in which 



charge transport is dominated by hopping processes. For 
the opposite case of a continuous-medium model, the dif- 
fusion coefficient is known to obey a generalized version 
of the Einstein relation; 25 ' 29 

where v(F) is the drift velocity. Our expressions for both 
RBM and REM, however, demonstrate that Eq. ([94]) is 
also violated in the general case, and the deviation is also 
proportional to We argue that, generally, there are 
no exact connections between D(F) and /J,(F) at F^O. 
Indeed, let us consider the REM. There are eleven quan- 
tities mi . . . mn (see Sec. IIV|) dependent on statistics 
of the energy levels and on the electric field. The mo- 
bility depends on three of them (7711,7713,7714), accord- 
ing to Eq. ([70]); the diffusion coefficient depends on all 
eleven quantities, see Eq. (|7ip . For an arbitrary den- 
sity of states, all the eleven quantities are independent of 
each other and cannot be reduced to each other. Conse- 
quently, there is no general way to reduce the diffusion 
coefficient to the mobility at non-zero electric field. 

In conclusion, we have examined analytically and 
numerically two models of one-dimensional hopping 
transport — the random-barrier (RBM) and the random- 
energy model (REM) . Exact analytical solutions of field- 
dependent diffusion coefficient have been obtained for 
both models in the case of nearest-neighbor hopping. 
We have demonstrated that the non-analytic field de- 
pendence ([9"Tj) of the diffusion coefficient, as well as the 
violation of the Einstein relation for any nonzero elec- 
tric field, are inherent properties of hopping transport in 
one-dimensional disordered systems. 
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